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Abstract 

We calculate the branching ratios and CP asymmetries of the B — )■ irp, ttlu decays in the per- 
turbative QCD factorization approach up to the next-to-leading-order contributions. We find that 
the next-to-leading-order contributions can interfere with the leading-order part constructively or 
' destructively for different decay modes. Our numerical results have a much better agreement with 

p, I current available data than previous leading-order calculations, e.g., the next-to- leading-order cor- 

r~| . rections enhance the — )• n^p^ branching ratios by a factor 2.5, which is helpful to narrow the 

gaps between theoretic predictions and experimental data. We also update the direct CP-violation 
parameters, the mixing-induced CP-violation parameters of these modes, which show a better 
agreement with experimental data than many of the other approaches. 
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I. INTRODUCTION 



The charmless B meson decays are not only suitable to study CP violations but also 
sensitive to new physics[l|. During the past decade, the B factory experiments achieved 
great successes. Furthermore, the current LHC experiments will provide 2-3 orders more B 
meson events than the B factories 2|. A large number of rare B meson decay channels will 
be measured by the future super B factories. The research on the charmless decays of B 
meson is therefore becoming more interesting than ever before 

The theoretical calculations of color-suppressed decay channels, such as B^ — )■ ir^ir^, met 
a difficulty for a relatively much smaller branching ratios than the experimental measure- 
ments The difference between direct CP-asymmetry measurement oi B^ — )■ K~^ii~ and 
B^ — 7- /i^7r° showed a very large discrepancy between the leading-order (LO) theoretical 
calculations and experimental data, which induced a lot of new physics discussions [3] . One 
of the standard model solutions to this puzzle also requires large color-suppressed tree am- 
plitudes Q|. Some of the next-to-leading-order (NLO) QCD calculations in the perturbative 
QCD factorization approach (pQCD) 8l-lll| show that the NLO contributions can signif- 
icantly change the LO predictions for some decay modes, especially the color-suppressed 
modes. It is therefore necessary to calculate the NLO corrections to those two-body charm- 
less B meson decays in order to improve the reliability of the theoretical predictions. 

The 5 — )■ vrp decays, which are helpful for the determination of the Cabibbo-Kobayashi- 
Maskawa(CKM) unitary triangle a angle measurement in addition to the i? — > tttt decays, 
have a much more complication. Either of B^ or B^ meson can decay to both the 7r^p+ 
and 7r+p~ final states, which lead to altogether four decay amplitudes. Since B^ and 5° 
meson mix easily, these channels exhibit unique features of mixing and decay interference 
in B physics. The recent B factory measurements indeed show that the interesting phe- 



nomenology with possible large direct CP asymmetry [12|. Unlike the branching ratios, the 
CP asymmetries are sensitive to high order contributions. Similar to the color-suppressed 
B^ — )■ 71^71^ mode, the neutral decay modes B^ — t- ir^p^, ir^u are also expected to receive 
considerable NLO contributions. Therefore, it is necessary to calculate NLO corrections to 
the B —J- vrp, vrcj decays in the pQCD approach for the reason that previous pQCD calcula- 



tions 



13| are already too old with only LO accuracy. In this paper, we calculate the NLO 



contributions arising from the vertex corrections, the quark loops and the chromo-magnetic 
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penguin operator O^g. Combining our results with the NLO accuracy Wilson coefficients 
and Sudakov suppression factors, we present a numerical analysis of -B — t- vrp, vro; decays. 

Our paper is organized as follows: we ffist review the pQCD factorization approach in 
Sec. mi Then, in Sec. IIIH we show our analytical results of NLO calculations. The numerical 
results are given in Sec. IIVI Finally we close this paper with a conclusion. 



II. THEORETICAL FRAMEWORK 



For the studied B — t- irp^TTUJ decays, the weak effective Hamiltonian He// for 6 — c? 
transition can be written as 

„ 10 

He// = -|[e„(Ci(/x)Or(/i) +C^2(/i)02"(/i)) 

^ j=3 

where C,u = VubV*^, = VtbV^*^ are the CKM matrix elements. Oi(/i) and C.t(/i) are the 
four-quark operators and corresponding Wilson coefficients, respectively. Expressions of Ci 
and Oi can be found in Ref.[l^. In the following, we will use this effective Hamiltonian to 
calculate decay amplitudes in the pQCD approach. So, we ffist give a brief review of pQCD 
approach and present relevant wave functions. 



A. pQCD factorization approach 

In the framework of the pQCD factorization, three scales are involved in the non-leptonic 
decays of B mesons: the weak interaction scale mw, the hard subprocess scale t and the 
transverse momenta of the constituent quark kj-. The large logs between W boson mass scale 
and the hard scale t have been resummed by the renormalization group equation method to 
give the effective Hamiltonian of four-quark operators. In two-body charmless hadronic B 
decays, the final state meson masses are negligible compared with the large B meson mass. 
Therefore the constituent quarks in the final state mesons are collinear objects in the rest 
frame of B meson. The momentum of light quark in B meson is at the order of Aqcd, such 
that a hard gluon is needed to transfer energy to make it a collinear quark into the final state 
meson. These perturbative calculations meet end-point singularity in dealing with the meson 
distribution amplitudes at the end-point. Usually in the collinear factorization approaches 
such as QCD factorization 4] and soft-coUinear effective theory [isl. people parameterize this 
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kind of decay amplitudes into free parameters to fit the data. While in the perturbative 
QCD factorization approach, we take back the parton transverse momentum kx to regulate 
this divergence. 

In the pQCD approach, the decay amplitude A{B — )■ M2M3) can be written conceptually 



as the convolution 



16| 



A{B ^ M2M3) = j d%d%d%Tr[CmB{ki)^M,{k2)^Ms{h)H{h,h,h,t)] (2) 

where ki are momenta of light quarks included in each meson, and Tr denotes the trace over 
Dirac and color indices. The hard function H{ki, k2, k^, t) describes the four-quark operator 
and the spectator quark connected by a hard gluon of order AMb, which can be calculated 
perturbatively. The energy scale t is chosen as the maximal virtuality of internal particles 
in a hard amplitude, in order to suppress higher order corrections ITj]. $1/- is the wave 
function of meson Mj. The hard kernel H depends on the processes considered, while the 
wave functions $a/. are process independent that can be extracted from other well measured 
processes, so one can make quantitative predictions here. 

It is convenient to work at the B meson rest frame and the light cone coordinate. The final 
state meson M2 is moving along the direction of v = (0, 1, Ot) and M3 is along n = (1, 0, Ot). 
Here we use Xi to denote the momentum fractions of anti-quarks in mesons, and kix to denote 
the transverse momenta of the anti-quarks. The mass of light meson (vr) is neglected. After 
integration over k^ , /c^ and k^ in Eq.([2]), we are led to 



A = J dxidx2dxzhidhih2dh2h^dh^ 

Tr[C(t)$B(xi, hi)^M,{x2, h2)^M,{x^, h)H{xi, bi, t)St{xi) exp(-5(t))] (3) 

where bi are the conjugate variables of kiT. The jet function St{xi) arises from the threshold 
resummation of the large double logarithms In^(xj), and the Sudakov exponent S{t) comes 
from the double logarithms of coUinear and soft divergences. 

B. Wave Functions 

There are generally two Lorentz structures in the B meson distribution amplitudes, which 



can be decomposed as [18 



With Nc = 3, they obey the following normalization conditions: 



^'V(fci) = -4=, /^0b(A.i) = O. (5) 



However, the contribution of (pB is numerically neglected [191] . Therefore, we will only con- 
sider the contributions from 0^. In b space the B meson wave function can be expressed 
by 

"^Bix, h) = -^={pB + mB)-f5(pB{x, h). (6) 

V ^^Vc 

For the light pseudo-scalar mesons tt, the wave function can be defined as [2^ 

$(P,x,0 = -^75[yP0^(x) +mo0^(x) +emoW " (7) 

V ^jVc 

where P is the momentum of the light meson vr, mo is the chiral mass which is defined using 
the meson mass irtp and the quark masses as mo = rn^p/imq^ + m^j). x is the momentum 
fraction of the quark (or anti-quark) inside the meson, respectively. When the momentum 
fraction of the quark (anti-quark) is set to be x, the parameter ^ should be chosen as +1(— 1). 

For the considered decays, the vector meson V{p^uj) is longitudinally polarized. The 
longitudinal polarized component of the wave function is defined as: 

= -^=y{mv(pvix) + pv^vix)) + mv(pv{x)], (8) 

where the polarization vector e satisfies Py ■ e = 0. 



III. ANALYTICAL CALCULATIONS 



Our NLO corrections for pQCD approach include the following parts: 

• The NLO hard kernel H^^\xi,hi,t), which includes the vertex corrections, the quark 
loops and chromo-magnetic penguins. 

• The NLO Wilson coefficients C^^^it), which have been calculated in the literature 

H. 

• The exponential Sudakov factor exp[— S'^'^*^(t)] includes the Sudakov factor s{P,b) 
and renormalization group running factor g2{t, b). 
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FIG. 1: NLO corrections to the hard kernels. The diagrams (a-f), (g,h) and (i,j) are commonly 
called vertex corrections, quark-loop corrections, and chromo-magnetic penguin corrections, re- 
spectively. 



So, at the NLO, Eq.Q can be written as 

A^^^ = j dxidx2dx3bidbib2db2b3db3Tr[C^^°{t)^B{xi,h)<^MAx2,b2)^Ms{x3,h) 

{H^'\x,, h, t) + H^'\x,, 6„ t))St{xi) exp(-5^^'^(t))]. (9) 



We will give these calculations in the following of this section. 



Vertex corrections 



The vertex corrections are part of the complete NLO Wilson coefficients for four-quark 
operators, which cancel the explicit renormalization scale n dependence of the Wilson co- 



efficients. The vertex correction diagrams are illustrated by 



'igsiT]|^a)-[Tl^f), among which 



Fig.(e) and (f) are new compared to the QCDF calculation Here, we have introduced 
transverse momentum kx in regularizing the infrared divergence. Our results are different 
from the QCDF approach ^ for different regularization schemes. 

The vertex corrections to the S — )■ Trp, iru decays modify the Wilson coefficients for the 
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emission amplitudes into 

ai(/i) ^ ai(/i) + '^Cpfj^Vr^M) + CMV[{M)l 

a.(/i) ^ a.(/.) + ^C^[^:^l-,(M) + a(/.)K(M)], ^ = 3 - 10, (10) 

where M denotes the meson emitted from the weak vertex, and the upper (lower) sign applies 
for odd (even) i. When the emitted meson M is a pseudo-scalar meson, the functions Vi{M) 
and V-{M) are given by 

l\n^-18 + ^-6t7T + ^J^dx<Pi{x)g,ix), for 1=1-4,9,10 
Vi{M) = { -16ln^ + 6 + ^ + ^-yff^J^dx^Ux)g2ix), for 1=5,7 

-16\n^ + 6 + ^ + ^^J^dx(j)^j{x)h{x) for 1=6,8 

-4 In ^ + ^ - Sivr, for 1=1-4,5,7,9,10 
y/(M) = ' ' (11) 

-16 In ^ + 12 + ^ for 1=6,8 

where ruh is the mass of b quark. The functions gi{x), g2{x) and h{x) are given as 

ln(l — x) In X 
^i(x) = 31na; + 31n(l -x) + ^ --2, 



X 1 — X 

+ [ln^x + 2Ll2( ) + 427rlnx - (x 1 - x)], 

X — 1 

ln(l — x) In X 
g^[x) = -31nx-31n(l -x) + 2- ^ ^ ' 



X 1 — X 

-Ffln^x + 2Ll2(— ^) -F427rlnx- (x 1 -x)], 
X — 1 

x Iri X 

h(x) = ln^x + 2Ll2( ) + —, r + 4z7rlnx - (x -> 1 - x). (12) 

X — 1 2(1 — x) 

When a vector meson V{V = p,co) is emitted from the weak vertex, (pMi'pAi) is replaced 
by 4>v{—(pv), and fm by fy in the third line of Eq.lfTT]). Note that, the amplitude 
from the operators O^s vanishes at LO, because neither the scalar nor the pseudo-scalar 
density gives contributions to the vector meson production, i.e. < + P\0 >= 0. On 
including the vertex corrections, the NLO piece aye, containing the vertex-correction of 05,8 
in Eq.f jTOj) . contributes through the following additional amplitudes jsl: 

fyF[^ ^ aycf^F[^ + fyF,^ (13) 
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where is the decay amphtude of factorizable emission diagrams with the structure of 
(y — A)^(y — A) insertion; while is the corresponding decay amphtude with (5 — P) ® 
{S + P) insertion. 



B. Quark loops 

The contributions from the quark loops are illustrated by FigJT|^g)-[^li). The quark-loop 
contributions are generally called the Bander-Silver-Soni mechanism 2]J], which plays a 
very important role in producing the direct CP-violation strong phase in the QCDF/SCET 
approaches. We include quark-loop amplitudes from the up-, charm-, and QCD-penguin- 
loop corrections, the quark loops from the electroweak penguins are neglected due to their 
smallness. 

For the h ^ d transition, the contributions from the various quark loops are described 
by the effective Hamiltonian "Hgjj 8], 

q=u,c q' ^ 

+ E %yt,V:,^^C^\^^^ /^)(J7p(1 - l,)T%){q'YTyi (M) 

q' ^ 

with 

= [GW(^,/2)-^]C3(/i)+ G^'^"\^,^')[c,{^^) + c,{^,)l (is) 

q" =u,d,s,c 

where P' being the invariant mass of the intermediate gluon, which connects the quark loops 
with the q'q pair. Because of the absence of the end-point singularities associated with 
/^,/'^ — 0, we have dropped the parton transverse momenta kx in /^,/'^ for simplicity. The 
integration function G^'^^fi, P) for the loop of the quarks q = {u, d, s, c) is defined as 

G(«)(/i,/2) = -4 / dxx{l-x)\n^ \ — . (16) 

Finally, the quark-loop contributions shown in Fig|T]^g) and [U^h) to the considered B — )■ 
ttV decays with V = can be written as 

^S.^ = {vm^:?f\B) = Yl Q-^vl + -^J^]- (17) 

q=u,c,t 



The two kinds of topological decay amplitude of the B V or B ^ n transition are written 

as 

= —i=Cpm'^^ / dxidx2dx3 / bidbib2db2(f)B{xi,bi) 
v3 Jo Jo 

X {[-(1 +X2)0y(x2)0;;^(x3) +2r^0y(x2)0f(x3) 

-(1 - 2x2)rv(Pt{xMUx2) + <PU^2)) 
+2(2 + X2)r^rv(j)v{x2)(l)^{x3) - 2x2r^rv(j)v{x2)(j)^{x3)] 
xal{h)hgi{xu X2, 6i, b2)C^'\h, f) exp[-S,i{h)] 
+2rv{2r^<P^{xs) - <P^{xs))<PU^2) 

y<al{t2)h,i{x2, xi, 62, bi)C^'\t2, /') exp[-S,i{t2)]} , (18) 

M'^Jv = -^C!%m% I dxidx2dxs / bidbib2db2(f)B{xi,bi) 
vS Jo Jo 

{ [(1 + X2)(p^ix2)<l)v{x3) + (1 - 2X2)r,(0^(x2) + 0nX2))0y (Xg) 
2ryr^(2 + X2)(j)^{x2)(j)v{x3) - 2ryX20^(x2)0y(x3) 

+2rv<P^{x2)<PU^3)] X a^(t^)/,^,(a;^, &2)C('')(^i, exp[-^,K^i)] 

+2r^[0f (X2)0^(X3) + 20f (X2)0t'(x3)] 

xa^s{t2)Ki{x2, xi, 62, &i)C('')(t2, /') exp[-S,i{t2)] } , (19) 

where the ratios ry = mv/mB,rn = rriQ/mB- The hard scales and the gluon invariant 
masses are given by 

ti = max{y/%mB, y/xiX2mB, \/ x^{l - X2)mB, 1/^2), 
t2 = max{y/x^mB, y/xiX2mB, ^/\x3 - Xi\mB, I/&1, 1/&2), 

= ^3(1 - X2)m|, /'^ = (x3 - xi)m|. (20) 

The hard functions hqi are included in the appendix. 
C. Chromo- magnetic penguins 

The chromo-magnetic penguin contributions are of NLO in within the pQCD formal- 
ism. They are at the same order in as the penguin contributions. 



According to ref. 22|, there are ten chromo-magnetic penguin diagrams contributing to the 



B decays, but only two of them are important, as illustrated by Fig. [Tt^i)and[T]|^j), while the 
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other eight diagrams are neghgible. The corresponding weak effective Hamiltonian contains 
the b ^ dg transition: 

= -^^^Cli-fOsg, (21) 

with 

Osg = ^m,d,a,.(l + 75)X^G"^''^&, (22) 
where i,j being the color indices of quarks. The corresponding effective Wilson coefficient 



The decay amplitudes of FigHJ^i) and[T|^j) can be written as 

M^^^^ = —=Cpm% / dxidx2dx3 / 6ic?6i62C?&2&3C?&3</'B(a^i, &i) 
v3 Jo Jo 

X {[(1 - X2)(Pt{x3)[2(j)v{x2) + ry(30^(x2) + 0U^2)) 

+X2rv{(pv{x2) - (Pv{x2))] - r^X3{l + X2)(30f (Xa) - 0^(x3))</)y(x2) 
-Vvil - X2){M^2) - 0U^2))(30r(^3) + 0^X3)) 
+VvX3{l - 2X2){^U^2) + (p'y{x2)){(pl{x3) - 30^(X3))] 

xC^^/a'^^{ti)hmgiA, B, C, bi, 62, b3)St{x2) exp[-S'„,g] 
+2ry[20;f(x3) + xsrMli^s) - 30^ (x3)]0t'(x2) 

xC'//a^,{t2)h„,g{A'^ B', C, 62, 61, b3)St{x,) exp[-5„3]} , (23) 

M.^^^ = —^Cpin^ / dxidx2dx3 / bidbib2db2b3db3(f)Bixi,bi) 
v3 Jo io 

X {[(1 - X2)0y(x3)[20^^(x2) + r^(30^(x2) + 0^(X2)) 
+a;2r,(0nx2) - (pv{^2))] + rvX3{l + X2)(30t.(x3) - (P'y{x3))(p^{x2) 

+Vv{l - X2){3MX3) + M^3)M{X2) - 0^X2)) 

+WX3(1 - 2x2)(30t^(x3) - 0Ua^3))(0^(a:2) + €{^2))\ 
xCf/al{h)h^g{A, B, C, bi, 62, b3)St{x2) exp[-Smg] 

+2r^[20y(a;3) - X3ry(0y(x3) - 30t^(x3)]0f (X2) 

xC|f a,2(t2)/im5(^', 5', C, 62, 61, 63)5t(xi) exp[-5„,] } , (24) 



where 



A = y/x^nif,, B = B' = ^/x^niB, C = ^/x^iY^^a^niB, 

A' = y/x^nib, C = - X3\mB. (25) 
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The hard scales ^1,^2 are the same as in Eq. fl20p . The hard function hmg and the Sudakov 
exponent Smg are given in the appendix. The jet function St{xi) can be found in Ref. |23|. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

Besides those specified in the text, the following input parameters will also be used in 
the numerical calculations jsl: 

iriB = 5.28GeV, r^o = 1.53ps, Tb+ = 1.638ps, 
fB = 0.21 ± O.OlGeV, \Vub\ = (3.47+°;}^) x 10-^ |Kd| = 0.97428, 
\Vtb\ = 0.999, \Vtd\ = (8.62;°;26) x 10-^ a = (90 ± 5)°. (26) 

The corresponding values of Aqcd are derived from asimz) = 0.1184 using LO and NLO 
formulas, respectively: 

LO: AgJ,j3 = (0-110 ±0.005)GeV, A^^l^ = O.USGeV; 
NLO : AgJ,j3 = (0.228 ± 0.008)GeV, AgJ^^ = 0.325GeV. (27) 

The B meson distribution amplitude is given by 

0b(x,6) = iVBx2(l -x)2exp[-— V - oM)']' (28) 
where the shape parameter Ub = 0.40 ± 0.04GeV has been fixed using the rich experimental 



data on the and i?^ decays p, uA 125 1 



For the vr meson, the twist-2 distribution amplitude (f)^{x), and the twist-3 distribution 
amplitudes 0^(x) and (j)^{x) are written as 13 1 



V 2iVc 

<P^{x) = ^^jl + 0A3Cl/\2x-l) + 0mCl/\2x-l)], 

<f^lix) = ^^(1 - 2x)[l + 0.55(10x2 - lOx + 1)] (29) 

with the pion decay constant /yr = 0.13GeV. The Gegenbauer polynomials are defined by 

Cl^\t) = ^(3t' - 1), Cy\t) = i(35t^ - 30t' + 3), 

Cl^'it) = li^t' - 1), C'J'it) = ^(21t^ - Ut' + 1) (30) 
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whose coefficients correspond to = 1.4GeV. 

The distribution amphtudes for the vector meson are hsted below 



13|: 



(f>vix) 



fT 

■'^ ■.[3{2x - If + 0.3(2x - l)2(5(2x - 1)^ - 3) 



2VQ 

+0.21(3 - 30(2x - ly + 35(2x - 1)^)], 

-^/^(l - 2x)[l + 0.76(10x2 - lOx + 1)], 
2v6 



(31) 



with the decay constant fp = 0.216GeV, fj = 0.165GeV, = 0.195GeV and 
O.MSGevfl. 



A. Branching Ratios 

The considered NLO contributions can interfere with the LO part constructively or de- 
structively for different decay modes. In Table [H we show our pQCD results for the CP- 
averaged branching ratios of the seven B ^ irp, nu decays together with the experimental 
data. In order to show the effects of the improvement, we use the same updated input 



paraments for the LO and NLO calcu 
than the previous pQCD calculations 



ations, which make the LO-pQCD predictions larger 



13|. Apparently, most of the NLO-pQCD predictions 



agree with the experimental measured values and better than the LO results. 



ac- 



For comparison, we also list theoretical predictions based on the traditional QCD 
torization approach (QCDF-I) 26[, modified QCD factorization approach (QCDF-II) 
which include the fitted penguin annihilation topology and color-suppressed tree amplitudes, 
and the ones obtained using SCET 28]. Comparing with the experimental data jol, it is 
easy to see that the LO-pQCD predictions are worse than the QCDF results, but our NLO- 
pQCD results have a better agreement with the experimental data. Our NLO predictions 
of the branching ratios for B — )■ vr^p^ decays are close to QCDF-II result but larger than 
those in SCET. Neglecting the small terms, it is due to the different B ^ tx and B ^ p 
form factors: SCET uses the smaller form factors F^'^'^ = 0.198 and Aq^^ = 0.291; while 
in our NLO calculations, F^^^ = 0.23 and Aq^^ = 0.30. 

For the color-suppressed tree dominant mode B^ — ?■ ir^p^, the NLO pQCD contribu- 
tions enhance its branching ratio by a factor 2.5, which are helpful to pin down the gap 
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TABLE I: Branching ratios (xlO ^) of i? — )• ttp,ttuj decays in the pQCD approach, together with 
results from the QCDF-J Q, QCDF-IlQ, 

the ones obtained from one solution of SCET and 

the experimental data 



Mode 


LO-pQCD 


NLO-pQCD 


QCDF-I [2^ 


QCDF-II [22] 


SCET [2^ 


Data 


B^B 




41.3 


r)c: 7+7.0+2.4+1.3+1.8 
'-5.5-1. 9-2.0-1. 6 


ofi [-+18.2+10.3+2.0+3.9 
OD.O_]^4 7_g g_3 5_2,g 


or 1+1.5+1.4 
^O.J-_2.2-1.8 


Ifi s+0.5+1.6 


23 ±2.3 


B+ ^ 




9.0 


r 4+1.4+0.5+0.6+0.3 
^■^-1.1-0.3-0.5-0.0 


11 Q+6.3+3.6+2.5+1.3 
^^•^-5.0-3.1-1.2-1.1 


o 7+2.7+1.7 
^■'-1.3-1.4 


7 Q+0.2+0.8 
'■^-0. 1-0.8 


8.3 ± 1.2 


B+ ^ 




14.1 


n f;+2.5+0.8+0.7+0.7 
-^■0-2.1-0.7-1.3-0.6 


14 n+6.5+5.1+1.0+0.8 
^^•^-5.5-4.3-0.6-0.7 


11 0+1.8+1.4 
11.0_1.1_1.4 


11 4+0.6+1.1 
^^■*-0.6-0.9 


10.9 ± 1.4 






0.15 


n O7+0.09+0. 02+0. 03+0.08 
^■'J' -0.08-0.01-0.05-0.02 


n /I+0.2+0. 2+0.9+0.5 
^•^-0.2-0. 1-0.3-0.3 


1 o+l. 7+1.2 
-^■•^-0.6-0.6 


1 c+O.l+O.l 
^■O-o. 1-0.1 


2.0 ±0.5 


B+ ^ 




8.4 


A fi+l. 2+0.5+0. 5+0.1 
^■"-0.9-0.4-0.4-0.1 


o 0+4.4+2. 6+1. 8+0.8 
''•''-3.5-2.2-0.9-0.9 


7+2.1+1.3 
"■'-1.0-1.1 


o C.+0.3+0.8 
'■°-0. 3-0.8 


6.9 ±0.5 






0.22 


n 00+0.06+0.01+0.04+0.04 
^■■^^-0.05-0.02-0.07-0.04 


n ni +0.00+0.02+0.02+0.03 
^•^^-0.00-0.00-0.00-0.00 


n ni +0.02+0.04 

"■"^-0.00-0.01 


n m i;+0.024+0.002 
U.UiO_o 000-0.002 


< 0.5 



between the pQCD calculations and the experimental data. This NLO BTZ{B^ ir^p^) is 
comparable with the result of QCDF-I, but still smaller than QCDF-II and SCET results 
and the experimental data. Soft corrections to a2 enhance the QCDF-II predictions, while 
in the SCET framework, the hard-scattering form factor (j is fitted to be relatively large 



and comparable with the soft form factor (. In a very recent paper [29[, the authors show 
the existence of residual infrared divergences caused by Glauber gluons in non-factorizable 
emission diagrams, which may resolve the large discrepancy between the theoretical predic- 
tions on BTZ{B^ — 7- TT^p^) and the data. For another color-suppressed tree dominant mode 
— i> Ti^uj, our pQCD prediction is comparable with the B^ Ti^p^ mode; while both 
QCDF and SCET predictions for this mode are less than B^ — )■ -K^p^ results. This should 
be clarified by future experiments. 

The theoretical uncertainties of the NLO-pQCD predictions are also shown in Table 
m The first error comes from the B meson wave function parameters Ub = 0.40 ± 0.04 
and Jb = 0.21 ± O.OlGeV; the second error arises from the uncertainties of the CKM 



matrix elements \V, 



ub 



(3.47: 



~0.16n 
-O.I2J 



X 10" 



td\ 



f8.62: 



h0.26\ 

-0.20 J 



X 10 ^ and the CKM angles 



a 



(90 ±5)°; the third error comes from the uncertainties of final state meson wave function 



parameters 



-0.1 II 
-0.2' "P 



0.18 ± 0.1 

,(5) 



20|; the fourth error is from the hard scale t 



aoogGeV, which characterizes the uncertainty 



varying from 0.75t to 1.25t and A^^,^ = 0.228 
of higher order contributions. It is easy to see that the most important uncertainty in 
our approach comes from the B meson wave function and CKM elements Vub- The total 
theoretical error is in general around 30% to 50% in size, which is smaller than the previous 
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leading-order calculation. 

Since both tree and penguin diagrams contribute to these decays, the decay amplitude 
for a given decay mode with b ^ d transition can be parameterized using CKM unitarity as 

A = CuT- = CuT[l + ze'^''+% (32) 

where the parameter z = \C,t/^u\\P/T\, the weak phase a = aTg[—C,t/C,u], and 6 = arg [P/T] is 
the relative strong phase between T and P part. The corresponding charge conjugate decay 
mode is then 

^ = - 6P = ^uT[l + ^e^(-"+^)]. (33) 

The CP-averaged branching ratio is 

Br{B^7Tp) = 1^1' + 1^1' ^ TB \^^T\^^i^2zcosacos5 + z\ (34) 

IGivmB 2 IGnmB 

which shows a clear CKM angle a dependence. This potentially gives a way to measure the 

CKM angle a by these decays, if we can really pin down the large theoretical uncertainties 

of the branching ratio calculations. For illustration, we show the LO and NLO results of 

Br{B^/B^ — i> TT^p^) in Fig [2] as a function of a with the hard scales varied from 0.75t to 

1.25t. We observe that the scale dependence of the NLO branching ratio is significantly 

smaller than that of the LO branching ratio, roughly from k, 50%, reduced to less than 10%. 



B. CP asymmetries 

Using fl32l) and fl33l) . we can derive the direct CP-violating parameter 

j^dir ^ ^ 22: sin a sin (5 

|^|2 + |^|2 l + 22cosacos5 + ^2- ^ ' 

It is clear that the non-zero direct CP asymmetry requires at least two comparable contri- 
butions with different strong phase and different weak phase. Since is proportional to 
sin a, it can be used to measure the CKM angle a, if we know the strong phase difference 
between the tree and penguin diagrams. The CKM angle a dependence of the direct CP- 
violating asymmetries of these decays are shown in Fig. [3l The accuracy of this measurement 
requires more precise theoretical calculation and more experimental data. 

The numerical results for the direct CP- violating asymmetries of B^ 7r^p°, p^vr^, vr^o; 
and B^ — )• 7r°p°, tx^oj decays are listed in Table [Tll The direct CP- violation parameters of 
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20 40 60 80 100 120 140 160 180 

a(degree) 



FIG. 2: The scale dependence of Br{B^/B^ t^^P^) of the LO(the black band) and the NLO(the 
gray band). 

—7- TT+p*^ is negative, while the direct CP-violation parameter of the other modes are 
positive. The direct CP-violation parameter of B'^ — t- tx'^uj is rather small for the almost 
canceled contributions of annihilation diagram, which are the dominant contributions to 
the strong phases in pQCD approach. Because the NLO Wilson evolution increases the 
penguin amplitudes and dilutes the tree amplitudes, the NLO direct CP- violation parameters 
(absolute value) of those decays are slightly enhanced compared with the LO predictions. 
However, for the color-suppressed tree dominant modes B^ — )■ 7r°p° and B^ — t- vt^cj, the 
direct CP asymmetry varies from —50% to 47% and from 52% to 98%, respectively. The 
big changes are attributed to a huge change of the strong phase of color-suppressed tree 
amplitudes caused by the vertex corrections. 

The theoretical uncertainties of the NLO-pQCD predictions are also shown in Table [III 
The first error shown in the table, comes from the B meson wave function parameters 
Ub = 0.40 ± 0.04 and J'b = 0.21 ± O.OlGeV; The second error arises from the uncertainties 
of the CKM matrix elements |Kfe| = (3.47+°;J^) x 10-^ \Vtd\ = (9.62+°;f) x lO'^ and the 
CKM angles a = (90 ± 5)°; the third error comes from the uncertainties of final state 
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TABLE II: The pQCD predictions for the direct CP-violating asymmetries of 



TT^p^, p^iT^,7r^(jj and B' 



in QCDF-I [261, QCDF-II 



decays(in units of %). We cite theoretical results evaluated 



27[, SCET 



28| and experimental data 



3 



for comparison. 



Mode 


LO 


NLO 


QCDF-1 [26j 


QCDF-Il [22 


SCET \28\ 


Data 






-26.4 


1 q 9+4.8+0.7+6.5+8.5 
^"•■^-5.3-0.7-5.5-9.6 


A 1+1.3+2.2+0.6+19.0 
^■-■■-0.9-2.0-0.7-18.8 


Q 0+3.4+11.4 
^•°-2.6-10.2 


1 Q 9+15.5+1.7 
^^■^-13.4-1.9 


18!?, 






20.1 


Y+4. 4+1.6+4.4+8.8 
-4.1-1.6-4.8-8.2 


4 r,+1.2+1.8+0.4+17.5 
^•"-1.2-2.2-0.4-17.7 


Q 7+2.1+8.0 
^•'-3.1-10.3 


19 n+9.4+0.9 
^'^•"'-lO.O-l.l 


2 ±11 






0.4 


r n+0.3+0.3+1.2+0.8 

'-'•"•-0.1-0.3-0.5-2.5 


1 O+0.5+2.7+0.8+2.1 
^•°-0.5-3.3-0.7-2.2 


iq 9+3.2+12.0 
io-'^_2.i_io.7 


9 q+13.4+0.2 
'^•■'-13.2-0.4 


-4 ±6 






-49.8 


4fi cr+8.4+2.2+6.2+7.1 
^"■"•-8.2-2.1-2.6-7.0 


ir; 7+4.8+12.3+11.0+19.8 
-■"••'-4.7-14.0-12.9-25.8 


11 n+5.0+23.5 
-■-■•'^-5.7-28.8 


o r+21.4+0.3 
''•"'-20.3-0.3 


-30 ± 40 






51.9 


07 f;+0.0+1.5+0.7+3.0 
■"-1.1-2.1-1.3-3.0 




17 n+55.4+98.6 
^'•"-22.8-82.3 


oq r;+79. 1+3.4 
■'^•"'-185.5-3.1 





0.18 ± 0.1; the fourth error is from 



meson wave function parameters 
the hard scale t varying from 0.75t to 1.25t and A^^^ = 0.228 -aoogGfeV, characterizing 
the uncertainty of higher order contributions. Unlike the CP-averaged branching ratios, the 
direct CP asymmetry is not sensitive to the wave function parameters and CKM factors, 
since these parameter dependence canceled out in Eq.(j35]). In addition, the CKM angles (a) 
uncertainty is quite small (~ 5%). Therefore, the most important uncertainties here are the 
scale dependence, which shows the importance of the NLO calculations. 



We also cite results evaluated in QCDF-I 



QCDF-II [273, SCET|28 



for comparison 



in Table [TTl Our predictions on direct CP asymmetries are typically larger in magnitude, 
most of which have the same sign with SCET approach. In QCDF framework, the strong 
phases are either at the order of or power suppressed in KqcD/f^b- So predictions in the 
QCDF-I approach on these channels are usually small in magnitude, most have different 



signs from our pQCD results 



271 quoted in Table HI 



In fact, the QCDF-II results 
already included large strong phase coming from penguin annihilation contributions, so that 
their results agree well with our pQCD ones. 

For the neutral decays, the situation is more complicated due to the B^ mixing. 

The CP asymmetry is time dependent 6| , when the final states are CP-eigenstates. A time 
dependent asymmetry can be defined by 

T{B\t)^f)-T{B\t)^f) 



T{B%t)^f) + T{B%t)^f) 
cos Ami, 



(36) 
(37) 
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a(degree) 



130 



FIG. 3: Direct CP-violation parameters of B — t- (the top band), — t- (the second 

band), S"*" — t- p^n^ (the third band), B^ — t- tt'^uj (the fourth band), B~^ — t- vr^'^p^ (the bottom 
band), function of CKM angle a 



where Am is the mass difference of the two mass eigenstates of the neutral B meson. The 
mixing-induced CP-asymmetry parameter Sf is referred to as 

2Imi\j) 



ir. 1 + ze 



i{&—a) 



(38) 

If penguin contribution is suppressed comparing with the tree contribution, we will have 
the approximate relation Sf ^ sin 2a for a negligible z parameter. From FigHJ one can see 
that the Sf is not a simple sin 2a behavior, since the z ~ 3.5 for n^p^ and z ~ 1.0 for tt^cu, 
reflecting a very large penguin contribution. 

The pQCD numerical results for the CP- violating parameters Sf of t^^P^, tt^w are 

displayed in Table UTTl together with the QCDF-II |27| and SCET 28|] results. It can be seen 
that the pQCD central value for S.,rQpO has a different sign from the other two approaches, 
because of the penguin contribution is bigger than the tree contribution in our approach. 
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TABLE III: The pQCD predictions for the CP-violating parameters Sj of B — )■ (in 



27| |. the ones obtained using SCET 



2m and 



unit of %), together with results from the QCDF-II 
the experimental data 6(]. The errors for these entries correspond to the uncertainties in the scale 
dependence and other input parameters, respectively. 



Mode 


LO 


NLO 


QCDF-II [27] 


SCET [28] 


Data [6] 


S^OpO 


47 


9^+26+9 
"^^-19-12 


9/1+15+20 
^^-14-22 


_iq+l4+lO 


10 ±40 




-37 


91+5+13 
^-•--10-11 


70+14+20 
'°-20-139 


79+36+7 
' ^-154-11 






FIG. 4: Mixing-induced CP-violation parameters S^OpO ( the gray band), mixing CP-violation 
parameters S^o^ (the black band), as a function of CKM angle a. 

Our theoretical errors for these entries shown in the table correspond to the uncertainties 
in the scale dependence and other input parameters, respectively. It is easy to see that the 
uncertainty is very large. Currently, no relevant experimental measurements for the CP- 
violating asymmetries of these decays are available. Our predictions for these quantities are 
different from those in QCDF-II and SCET. We have to wait for the experimental data to 
resolve these disagreements. 
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TABLE IV: The LO-and NLO-pQCD predictions for the CP-violating parameters C, S, AC and 
AS* of B^/B^ — > TT^pf (in units of^%), together with results from the QCDF-I ' ^ ^ 

and the experimental data ^. 



the ones obtained using SCET 



261, QCDF-Il|27l, 



The errors for these entries 



correspond to the uncertainties in the scale dependence and other input parameters, respectively. 



Mode 


LO 


NLO 


QCDF-I [26j 


QCDF-II [27] 


SCET [28] 


Data [6] 


Acp 


-11 


, 7+4+4 
-■^'-3-4 


1 +0+1+0+10 
J^-0-1-0-10 


_ii+0+7 
-'--'^-0-5 


91+3+2 
^-•■-2-3 


-13 ±4 


C 


6 


1 c;+2+2 


n+0+1+0+2 

'-'-0-1-0-2 


Q+0+5 
^-0-7 


1+9+0 
^-10-0 


1 ± 14 


s 


-12 


01+6+16 
"^^-3-15 


1 q+60+4+2+2 
^"^-65-3-1-1 


A+l + W 

^-1-9 


1+6+8 
^-7-14 


1±9 


AC 


17 




1f;+6+23+l + l 

iU_7_26-2-2 




1 9+9+1 
^^-10-1 


37 ±8 


AS 


-7 


7+0+2 
'-0-1 


9+1+0+0+1 
^-0-1-0-1 


9+0+3 
^-0-2 


40+5+3 


-5± 10 



C. Time dependent asymmetry parameters of B^{B^) — )• vr^p^ decays 



Both and can decay into both the ir^p^ and final states. This is an interesting 
example of CP asymmetry in B decays, which is the only measured combination of four 
channels. Aj, Af, Aj and Aj are defined as follows 3o| : 



Af = (7r-p+l^e//|5°), Af-={7r+p-\H,ff\B'), 

Af = (7rVl^e//|5°), A;=(7rVl^e//|5°). 



(39) 



The system of four decay modes can define the time- and fiavor-integrated charge asymmetry: 



A 



CP 



\Af 


2 + 


Af 


2 _ 




2 _ 




2 


\Af 


2 + 


Af 


2 + 




2 + 


Af 


2 



(40) 



In the standard approximation, which neglects CP violation in the B^ — B^ mixing matrix 
and the width difference of the two mass eigenstates, the four time dependent widths are 
given by the following formulas 13| : 



r(5°(t) ^ 7r-p+) = e-^'^{\Af\'^ + \Af + Cf cos Amt- Sf sin Amt], 



r(5°(t)^7rV) = e~^'-i\Aj-\ 



\Af\^)[l - Cf cos Amt + Sf sin Amt], 



-rt 



T{B%t) ^ n+p-) = e~^'^{\Aff + \Aff )[! + Cf cos Amt- Sf sin Amt], 



e ^"2(1^/! + l^/r)[l - Cf cos Amt + Sf sin Amt], 



(41) 
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where Am > denotes the mass difference, and T is the common total width of the B meson 
eigenstates. Cf and Sf are defined as 

_ \Af\'-\Af\' _ 2Im{\f) _ 6 Af 

- w+\Afr ^'-i + \Af/Afr ^''^ 

For decays to the CP-conjugate final state, one replaces / by / to obtain the formula for Cj 
and Sf. Furthermore, we define C = |(C/ + Cj), S = j{Sf + Sf), AC = ^{Cf — Cj) and 
AS* = \{Sf — Sf). S is referred to as mixing-induced CP asymmetry and C is the direct 
CP asymmetry, while AC and AS* are CP-even under CP transformation Aj— j-l/Aj. If/ 
is CP eigenstate there are only two different amplitudes since / = /, and AC, /S.S vanish. 
The complicated formulas (jH]) return back to the simpler one in Eq. fl37|) . 
According to ( j32l) and ( j33l) . we can write Eq.( j39l) as 

Af = iuT-itP. A; = e.T'-6P', 

Af = CuT' - QP', Af = CT - ^P, (43) 

where T and P denote the tree diagram amplitude and penguin diagram amplitude of — >■ 
p^ii~ , respectively; while T' and P' denote the tree diagram amplitude and penguin diagram 
amplitude of — t- 7r+p", respectively. The asymmetries AS* ~ |y|2^|y/|2 sin 6 cos a are 
suppressed by the small penguin-to-tree ratios (|P/T|, |P'/T'| ^ 1) and the small relative 
phase 9 between T and T' {9 ^ 3.4°), hence they are always small in pQCD factorization. 



This conclusion is similar to that in QCDF 



26 



27|, although the absolute magnitude of A^* 



are much larger in pQCD 
decays including the LO 



than in QCDF. All the CP- violation parameters of B^ / B^ ix^p^ 



QCDF-II |27l, SCET 



13| and NLO results of pQCD, QCDF-I 
and the experimental data are collected in Table IIVI It is clear that the NLO-pQCD 
prediction for the CP- violation parameter Acp, AC and AS* agrees with the experimental 
results very well. The predictions of pQCD for CP-violation parameters in Table |IV] are 
comparable with the QCDF-II, and are better than QCDF-I and SCET predictions, which 



is also shown in other B decay channels 



V. CONCLUSION 



|3l|. 



In the framework of the pQCD approach, we calculated the NLO QCD corrections to 
the B — 7- vrp, txoj decays including the vertex corrections, the quark loops, the magnetic 
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penguin, and the NLO Wilson coefficients, the Sudakov factor and RG factor. We found 
that the NLO corrections improved the scale dependence significantly, and had great effects 
on some of the decay channels. Our NLO-pQCD calculations agree well with the measured 
values. For example, compared with LO predictions, the NLO corrections decease (increase) 
the branching ratio of B^/B^ — t- 7r^p^(i?° — > tt^p^), and improve the consistency of the 
pQCD predictions. The NLO corrections play an important role in modifying direct CP 
asymmetries. For the color-allowed tree dominant modes, the NLO Wilson coefficients 
enhance the penguin amplitudes, the larger subdominant penguin amplitudes increase the 
magnitudes of the direct CP asymmetries due to the stronger interference with the dominant 
tree amplitudes. The predictions of pQCD for CP-violation parameters are better than 
QCDF-I and SCET predictions. 
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We show here the hard function hgi and h^g the Sudakov exponents Sqi^mg{t) appearing 
in the expressions of the decay amplitudes in IIIIl 



Appendix 



hqi{xi,X2,bi,b2) = Ko{y/xiX2mBbi) 



X [9{b 



'1 



' b2)KQ{^Jx^mBbi)lQ{y/x^mBb2) 
bi)Ko{^/x^TnBb2)Io{y/x^mBbi)]St{x2), 



+9{b2 



(44) 



hmg{A,B,C,bi,b2,b3 




Jo (AbitanO) Jq {Ab2 tan 9) Jo {Ab-s tan 9) 



(45) 



21 



where Jo is the Bessel function and Kq, Iq are modified Bessel functions with Ko{—ix) = 
-{ti/2)Yo{x) + i{tx/2)Jo{x). 

The Sudakov exponents used in the text are defined by 

5 

Sqi{t) = s{ximB, bi) + s{x2mB, ^2) + s{{l - X2)mB, ^2) + ^92{t, bi) + 2g2{t, 62), (46) 

Smgit) = s{ximBM) + s{x2mB,b2) + - X2)mBM) + s{x^mBM) 

- X3)mB, 63) + \92{t, 61) + 2(?2(t, b2) + 2(72(t, ^3) (47) 



where the functions s(P, b) have been defined in Ref. 32|. The RG factor (72 (^, b) is given by 



2 ln(t/AgcD) . , 2/3, l^a^(pAfci;)) ^"^^"^Afc)) 
Po'^'^-HbAccoV^ In(^) 

1 1 

+ 
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/3o = 11 -^'^Z, /3i = 102-yn^ (48) 
where nj is the number of quarks with mass less than the energy scale t. 
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